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ABSTRACT 

Within the cosmic framework clusters of galaxies are relatively young objects. Many 
of them have recently experienced major mergers. Here we investigate an equal mass 
merging event at z w 0.6 resulting in a dark matter haloe of ~ 2.2 x 10 14 /i _1 M at 
z = 0. The merging process is covered by 270 outputs of a high resolution cosmological 
N-body simulation performed with the ART (adaptive refinement tree) code. Some 
2 Gyrs elapse between the first peri-centre passage of the progenitor cores and their 
final coalescence. During that phase the cores experience six peri-centre passages with 
minimal distances declining from ~ 30 to ~ 2 h~ kpc. The time intervals between the 
peri-centre passages continuously decrease from 9 to 1 x 10 8 yrs. We follow the mean 
density, the velocity dispersion and the entropy of the two progenitors within a set of 
fixed proper radii (25, 50, 100, 250, 500, 1000 /i _1 kpc). During the peri-centre passages 
we find sharp peaks of the mean densities within these radii, which exceed the sum of 
the corresponding progenitor densities. In addition to the intermixing of the merging 
haloes, the densities increase due to contraction caused by the momentary deepening 
of the potential well. The velocity dispersions also peak during peri-centre passages. 
Within the fixed proper radii the entropy of the most massive progenitor after the 
merger settles close to its pre-merger values. At the end of the oscillatory relaxation 
phase the material originating from the less concentrated of the two equal mass pro- 
genitors is deposited at larger radii and shows a slightly more radially anisotropic 
velocity dispersion compared to the material coming from the more concentrated pro- 
genitor. Every peri-centre passage is accompanied by a substantial drop of the central 
potential well. We briefly discuss the possibility that AGN outbursts are triggered by 
the periodically changing potential. 
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1 INTRODUCTION 

The current cosmological model (e.g. ISpergel et al.l 120061 ) 
assumes the gravitating matter to be composed of ~ 15 
per cent baryonic and ~ 85 per cent dark matter (DM). 
The former interacts electromagnetically, the latter is not 
directly visible. By comparing lensing observat ions with X- 
ray m aps of the merging cluster 1E0657-558 IClowe et al.l 
(|2006t ) find strong evidence for unseen matter concentra- 
tions that are more massive than the baryonic component. 
These observations emphasise the complexity of the interac- 
tions between dark and baryonic matter during mayor merg- 
ing events. Based on hydrodynamic si mulations of head-on 
mergers between model galaxy clusters IPearce et al.1 (|1994T ) 
find a damping of the core oscillations by the gaseous com- 
ponent. The situation may become even more complex if 
the collision of gas-rich systems can trigger star formation 
and outflows from the central acti ve galactic nucleus (AG N). 
Multi wave band observations bv lConselice et all (|200lh re- 



veal that this may actually the case for the central galaxy 
NGC 1275 in the Perseus Cluster. 

Using high resolution gas dynamical simulations 
iNagai et al.l (|2006l ) demonstrate that the X-ray observations 
of unrelaxed clusters cannot accurately measure the density 
and temperature gradients in the hot baryonic intra cluster 
medium (ICM). In addition some potentially relevant phys- 
ical processes, such as AGN bubbles, magnetic fields, and 
cosmic rays, are not well understood yet and are usually ex- 
cluded from such simulations. Here we do not aim to model 
the behaviour of the gas directly. Our analysis is based on 
a pure N-body simulation. However, since the collisionless 
component dominates the gravitational potential, we hope 
to gain some complementary insight into the violent pro- 
cesses during a mayor merger of galaxy clusters. The merg- 
ing event is covered by a set of 270 simulation outputs which 
allows us to precisely measure even highly variable features 
during the violent relaxation phase. Our goal is twofold: (1) 
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Table 1. Most massive FoF groups for linking lengths of 11 = 
0.20, 0.17, 0.10, 0.07, 0.02 times the mean particle distance at a 
redshift of z = 1.2. The three columns indicate the membership 
of the (sub)structures to three spatially separated haloes. The 
omitted entry for the 0.20 groups means that the largest 0.20 
group includes the two largest 0.17 groups below. 

we provide a precise time table of events during the core 
merging phase. The dynamics may be slightly altered by 
the presence of collisional material, however the general be- 
haviour is not expected to change profoundly. (2) we in- 
vestigate the evolution of the central potential to achieve 
some basic understanding of the gravitational forces which 
the baryons are exposed to. Ultimately we wish to uncover 
interactions between the DM and the baryonic components 
which result in observable phenomena. 

The paper is organised as follows. § [2] introduces the 
simulation and data preparation. In §[3] we present the evo- 
lution of the density, the velocity dispersion and the entropy 
of the DM within spheres of fixed proper radii. Additionally, 
the evolution of the central gravitational potential is dis- 
cussed. Finally, §[4] provides a summary of the main results 
and a discussion focusing on the potential interplay between 
the dark and the baryonic component. 



2 SIMULATION AND DATA PREPARATION 

The N-body simulation of a cluster-sized DM haloe has 
been performed with the a daptive mesh refinement code 
ART (jKravtsov et al.lll997h . The present analysis is based 
on ~ 270 snapshots between the redshifts z = 1.2 and z = 0. 
The cluster haloe has been selected from an initial low- 
resolution simulation, 128 3 particles within an 80/i -1 Mpc 
cube in a ACDM universe ( fi m = 0.3; D.a = 0.7; h = 0.7; 
as — 0.9). Higher mass resolution ha s been achieved by 
means of the multiple-mass technique l|Klvpin et al] 120011 ) 
applied to the particles within a Lagrangian volume of ap- 
proximately twice the virial radius (2 x 1.23 /i -1 Mpc) of the 
cluster at z = 0. This technique results in a mass resolution 
of 3.2 x 10 s /i _1 Mq which is equivalent to an effective resolu- 
tion of 512 3 particles within the 80/i _1 Mpc box. The mini- 
mum cell size allowed in the ART runs was set to 1.2 ft _1 kpc. 
The simulation started at z — 50 and the virial mass of the 
cluster at z — has is ~ 2.2 x 10 14 /i _1 Mq corresponding to 
~ 700000 particles within the virial radius. The haloe anal- 
ys ed here is the most mas sive cluster in the sample discussed 
bv lAscasibar et all (|2004h . 

At z ~ 0.6 two progenitors of similar mass merge. They 
are by far the largest units assembled by the cluster up to 
the present time. We wish to identify those particles at later 
times which belonged to one of the two most massive pro- 
genitors before the merger. To that purpose we apply a hi- 
erarchical friends-of-friends approach (FoF, see lKlypin et al] 
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Figure 1. Density profiles of the merged cluster at z = 0.2 and 
the two most massive progenitors at z = 0.85. The inner verti- 
cal lines display the characteristic radii of the NFW-fits and the 
outer vertical lines show the virial radius. The humps seen in both 
progenitor profiles at ~ 700 /i~ 1 kpc indicate the locations of the 
other progenitor at that redshift (compare to Fig. [3J. 

Il999l ) at z = 1.2, which is ~ 2 Gyr before the actual merger 
takes place. First the three most massive haloes within the 
high resolution Lagrangian volume are found by a linking 
length of II — 0.17 times the mean particle separation. With 
a sequence of decreasing linking lengths, namely 0.10, 0.07 
and 0.02 times the mean particle separation, substructures 
of increasing density are detected. Tab.[T]lists the three most 
massive structures for the various linking lengths. 

We find that the three most massive groups detected 
for each linking length are associated one by one with the 
three most massive spatially separated 0.17 haloes. Visual 
inspection reveals that the second most massive progenitor 
is in the final state of a merging event at 2 = 1.2. However, 
at z fa 0.85 when the central parts of the two progenitors 
begin to penetrate each other, the relaxation associated with 
that precedent merging event is entirely completed. There- 
fore, the two progenitors can be assumed to be fairly undis- 
turbed objects at the onset of the major merger. In the fol- 
lowing we will refer to the two most massive progenitors, 
comprising 153916 and 140898 particles, as the first and 
second ranked progenitors. The difference in mass is only 
~ 9 per cent, which justifies the perception of an equal mass 
merging event. To get an idea of the constellation relative to 
each other we apply a further FoF analysis this time with an 
enlarged linking length of 0.20 times the mean particle sepa- 
ration listed in the first line in Tab.[T] The most massive 0.20 
group embeds the two most massive 0.17 haloes, indicating 
that those two are about to merge. The second 0.20 group 
is much less massive and encloses the third 0.17 group, indi- 
cating that the third 0.17 group is well separated from the 
two others. In fact this group will merge with the main clus- 
ter at z ~ 0.1, a long time after the completion of the major 
merger. To summarise, at z = 0.85 the two most massive 
progenitors are fairly relaxed systems. The second ranked 
progenitor shows disturbances due to a precedent merging 
event at z > 1.2, however these are completely erased at 
2 = 0.85. The two progenitors have similar masses and are 
about to merge. The third most massive object is clearly 
separated at 2 = 1.2 and merges with the main cluster not 
before 2 = 0.1. 

Fig. [T] compares the density profiles of the first and 
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Figure 2. The orbits of the progenitor cores projected onto the 
xz-plane (left panel) and the yz-plane (right panel). The orbits 
are plotted between z ~ 0.9 and z ~ 0.3. 
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Figure 3. Evolution of the distance between the two progenitor 
cores with time. The vertical lines indicate the 6 succeeding peri- 
centre passages. 



second ranked progenitors at z — 0.85 with the density 
profile of the merged system at z — 0.2. In agreement 
with the FoF analysis, which indicated that the second 
ranked progenitor formed later, the first ranked progenitor 
shows higher concentration (c = 6.99) than the second one 
(c = 6.57), where the concentration results from a NFW-fit 
|Navarro et al.lll997l ) to the density profiles. The densities 
within the merged system are about two times larger than 
the densities within the two most massive progenitors. The 
hump seen in both progenitor profiles at r ~ 700 h~ kpc 
reflects the location of the other progenitor. 

Simulations provide an opportunity to trace the evolu- 
tion of single particles. This is of particular interest since it 
enables us to analyse the properties of the subset of parti- 
cles within the merged cluster that originate from each pro- 
genitor. In addition the orbits of the progenitor remnants 
within the already merged system can be easily traced if the 
substructure identification is based only on those particles 
which originate from the corresponding progenitor. There- 
fore, we keep record of the particles which belonged to one 
of two most massive haloes at z = 1.2. That way the pro- 
genitor of a given particle can be identified for any given 
subsequent time. A particle can either originate from the 
most or the second massive progenitor at z — 1.2 or from 
neither of them. 

To trace the orbit of a progenitor during the merging 
process we compute the top hat smoothed density field, 
Ttophat = 25/i _1 kpc (proper), based only on those parti- 
cles which belonged to the corresponding 0.17 FoF group 
at z = 1.2. Subsequently, we determine the location of the 
peak density. Applied to all snapshots, this procedure en- 
ables us to trace the orbit of the progenitor core during 
peri-centre passages and at late evolutionary states even if 
there is a strong contamination of particles from the other 
progenitor. Fig. [2] displays the projected orbits of the two 
most massive progenitors between z ~ 0.9 and z ~ 0.3. A 
complementary aspect is provided by Fig. [3] which shows the 
evolution of the distance between the two progenitor cores 
with time. During the first peri-centre passage the minimum 
distance of the progenitor cores is ~ 30/i _1 kpc and it grad- 
ually declines to ~ 2 ft _1 kpc for the last resolved peri-centre 
pa ssage. Similar d e cays o f su bhaloe orbits have been show n 
bv lKravtsov et al.l (|2004) and lBovlan-Kolchin fc Mai (|2006l) . 
Our approach enables to trace six peri-centre passages before 
the two progenitor cores entirely fuse to one single structure. 



About 2 Gyr elapse between the first peri-centre passage and 
the final coalescence. 



3 RESULTS 

The present analysis of an N-body merging event takes ad- 
vantage of the possibility to mark particles according to their 
association with one of the two most massive progenitors 
at z = 1.2. Consequently the properties of all the parti- 
cles within the final halo volume can be compared to the 
properties of only those particles that originally were asso- 
ciated with either the first or the second ranked progenitor. 
In the foregoing section we have verified that the masses 
of the two progenitors at the onset of the merging process 
(z ~ 0.85) are in fact very similar but the concentrations 
differ. The first ranked progenitor is somewhat more con- 
centrated. This property is important for the interpretation 
of the results presented here. The following paragraphs fo- 
cus on the evolution of mean values of the density, velocity 
dispersion and the entropy within a variety of proper halo- 
centric radii (r = 25, 50, 100, 250, 500 and 1000 /T 1 kpc). 
Finally we discuss the oscillations of the central potential 
during the merging process. 



3.1 Density 

The evolution of local densities during a merging event of 
collisionless systems depends on two competing processes: 
(1) the volume overlap of the two systems causes an increase 
of the local densities; (2) the deepening of the gravitational 
potential well that induces a contraction accompanied by 
accretion of additional material from the surrounding. To 
disentangle these two mechanisms Fig. [4] shows the evolu- 
tion of the mean densities in two different ways. The left 
hand side of Fig. [4] displays the evolution of the 'total' den- 
sities, i.e. densities based on all particles within the corre- 
sponding halo-centric spheres independent of the particles' 
origin (first or second ranked progenitor or none of them). 
The right hand side of Fig. [4] exhibits the densities based 
only on those particles which originally (at z — 1.2) were 
associated with each progenitor. In both figures the heavy 
lines are used for the density evolution associated with the 
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Figure 4. Left: Evolution of the mean densities for the two progenitors within the listed radii scaled by the mean densities of the first 
ranked progenitor (pi) at z = 0.2. The horizontal dotted lines mark a value of 1 for the pair of graphs next to it and a value of for the 
pair of graphs above. The heavy lines trace the evolution of the first ranked, more concentrated progenitor. The thin lines display the 
evolution of the second ranked, less concentrated progenitor. The vertical lines indicate the 6 succeeding peri-centre passages as introduced 
in Fig. [3] Right: Evolution of the mean densities for the two progenitors within the listed radii based only on the particles which were 
associated with the related progenitor at z = 1.2. The densities are scaled by the mean densities of the first ranked progenitor (pi) at 
z = 0.2. Heavy lines indicate the evolution of the first ranked, more concentrated progenitor and thin display the behaviour of the second 
ranked, less concentrated progenitor. The vertical lines indicate the 6 succeeding peri-centre passages as introduced in Fig. [3] 



first ranked progenitor and the thin lines depict the evolu- 
tion of the second ranked progenitor. All densities are scaled 
by the mean densities of the first ranked progenitor (pi) at 
2 = 0.2. (In fact, if the densities are computed based on all 
particles within a sphere (pi) = {p^j since two progenitor 
cores coincide at z — 0.2). Each horizontal dotted line indi- 
cates a value of one for the pair of graphs next to it and a 
value of zero for the pair of graphs above. 

For all proper inclusion radii the evolution of the mean 
densities in the left panel of Fig. [4] passes a maximum value 
around z = 0.7. For all radii lesssim500 /i -1 kpc the mean 
densities after the completion of the merging are somewhat 
lower than the sum of two the corresponding progenitor den- 
sities. For example the tracks of the mean densities within 
100/i~ 1 kpc raise from ~ 0.6 for both progenitors to ~ 1.5 
at the first peri-centre passage. After passing a subsequent 
minimum and a less pronounced second maximum the den- 
sity of the merged system settles at a value of 1, whereas the 
sum of the two initial densities results in a value of 1.2 . In 



contrast to the final density, the density at first peri-centre 
passage (~ 1.5) is larger than the sum of the two initial den- 
sities. This additional increase cannot be explained solely 
by the volume overlap of the two progenitors, it must be 
at least partly caused by the deepened potential well during 
the peri-centre passage which attracts an additional amount 
of matter towards the centres of the progenitors. 

The right hand side of Fig. U measures the density en- 
hancement as a response of the deepened potential well dur- 
ing peri-centre passages. Here the computation of the den- 
sities is based only on the particles which were associated 
with each corresponding progenitor at z = 1.2. This ap- 
proach rules out increasing densities as a result of overlap 
with the other progenitor. The only way that the mean den- 
sities within the fixed proper radii can increase is by contrac- 
tion, causing additional material to stream inwards across 
the boundary of that particular sphere. The densities are 
normalised by the mean densities associated with the first 
ranked progenitor (pi) at z = 0.2. In all cases the mean 
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densities are lowered after the merging event and the den- 
sities of the second ranked and less concentrated progenitor 
decrease by somewhat higher fractions. Some of the kinetic 
bulk energy of the progenitors is used to accelerate individ- 
ual particles which subsequently leave the original volume. 
The modest decrease for the 1000/i -1 kpc is explained by 
the fact that before merging the two most massive progen- 
itors have virial radii ~ 500ft _1 kpc (see Fig. [TJ, similarly 
the extend of the 0.17 FoF groups at z = 1.2 is substan- 
tially smaller than the 1000/i _1 kpc spheres. Therefore, if 
particles are pushed towards larger radii during the merging 
event they may still reside within the 1000/i _1 kpc sphere. 
However even for the 1000 /i kpc spheres a 10 percent loss 
can be observed. 

It is of interest to see how much additional matter (i.e. 
particles) not associated with either of the two most massive 
progenitors at z — 1.2, can be found within the various 
spheres after the merging event. In the following table we 
display the sum of the particles descending from both most 
massive progenitors divided by the total number of particles 
within the listed radii. 

h^kpc 25 50 100 250 500 1000 

(Ni+N 2 )/N 1.00 0.99 0.95 0.84 0.68 0.52 

The table above indicates that all particles within the central 
25 /i _1 kpc and 95 per cent of the particles within 100 h~ kpc 
descend from one of the two most massive 0.17 FoF groups 
at z — 1.2. Within the 1000 ft - kpc sphere ~ 50 per cent 
of the material was not associated with either of the two 
progenitors. As discussed above this is mainly caused by the 
fact that the 1000 h~ kpc sphere is substantially larger than 
the FoF groups, therefore already at z = 1.2 a substantial 
fraction of particles is included which is not associated with 
either of the two 0.17 FoF groups. 

The vertical lines in Fig. U indicate the moments of 
peri-centre passage derived from the evolution of the core 
distances in Fig. The right hand side of Fig. U reveals a 
slightly delayed peak of the densities based only on the orig- 
inal progenitor particles compared to the densities based on 
all particles within the given radii. The density enhancement 
due to the overlap of the two progenitors is followed by a 
density enhancement due to a contraction of the original par- 
ticles. For example the first peaks for the 50 ft _1 kpc sphere 
on the right hand side of Fig. 2] are delayed by ~ 100 Myr 
compared to the corresponding peaks on the left hand side. 
For larger inclusion radii this delay time increases. As clar- 
ified in Fig. [3]the 25 h~ 1 kpc spheres of the two progenitors 
overlap only partially during the first peri-centre passage 
with a minimum core distance of ~ 30/i _1 kpc. Thus most 
of the density increase within the innermost sphere must 
be caused by contraction involving mostly original parti- 
cles, which should erase differences between the two panels. 
Despite of the slightly higher density peak for the second 
ranked progenitor (left panel) and a marginal delay for both 
peaks of the individual progenitors (right panel), the ampli- 
tude and shape of the first peaks of the 25/i _1 kpc spheres 
are very similar in both panels. 

The evolution of the densities has revealed that the peak 
densities during peri-centre passages within fixed proper 
radii are caused by the volume overlap and additionally by 
contraction as a response of the deepened potential well dur- 



ing that phase. The final masses included within the fixed 
radii are not the sum of the masses within the corresponding 
radii of the progenitors. The transformation of kinetic bulk 
energy into random velocities pushes individual particles out 
of the original volume. Despite having very similar masses, 
the lower concentrated progenitor loses a higher fraction of 
its original central particles compared to the higher concen- 
trated progenitor. 

3.2 Velocity dispersion 

Using an approach similar to that of the densities also the 
mean velocity dispersions are investigated in a twofold man- 
ner. Fig [5] shows the evolution of the mean velocity disper- 
sion for all particles within the corresponding spheres and 
for only those particles which were associated with each in- 
dividual progenitor at z — 1.2. The heavy lines display the 
evolution of the first ranked progenitor whereas the thin 
lines display the fate of the second ranked progenitor. Be- 
fore computing the velocity dispersions the bulk velocity has 
to be subtracted from the individual particle velocities. For 
the the bulk velocity of the merging subhaloes we adopt the 
mean velocity within the central 25/i _1 kpc sphere, consid- 
ering only those particles associated with each progenitor at 
z = 1.2. This approach avoids a high background contami- 
nation, in particular during peri-centre passages. In addition 
the small volume sphere (r = 25/i _1 kpc) allows us to focus 
on the very central processes. As a test we also computed the 
evolution of the velocity dispersions based on bulk velocities 
within 250 h. -1 kpc spheres resulting in negligible deviations. 
The red vertical lines depict the moments of the peri-centre 
passages obtained from Fig. [3] 

We find the velocity dispersions within fixed proper 
radii do not change by more than 20 per cent between the 
values before and after the merging event. This is valid for 
the total velocity dispersion displayed in the left panel of 
Fig. as well as for the velocity dispersions based only on 
individual progenitor particles on the right. During the first 
peri-centre passage the velocity dispersions exhibit a mo- 
mentary increase of < 60 per cent. At the subsequent apo 
centre passage the velocity dispersions recede to the pre- 
merging values. The dispersions within the 25 /i _1 kpc sphere 
drop even slightly below the initial values. At the following 
peri-centre passage a second peak with substantial lower am- 
plitude is reached. The following oscillations for the total 
velocity dispersions indicate some relaxation process but a 
synchronisation with the peri-centre passages is ambiguous. 
A slightly different picture arises for dispersions based only 
on the original progenitor particles (right panel of Fig. [5|. 
At least for the central dispersions (r < 100ft _1 kpc) the 
particles descended from the second ranked progenitor show 
a synchronisation with the peri-centre passages. 

The double peak structure in both panels of Fig. [5] in 
particular for the spheres with intermediate radii resem- 
bles the behaviour of the total kinetic energy found by 
IPearce et all (| 19931 ) who analysed isolated N-body merg- 
ers. However, one has to keep in mind that our approach 
probes the evolution of the mean velocity dispersion along 
the paths of the progenitors which is not exactly the same 
as the total kinetic energy of the whole system. Neverthe- 
less, the fact that we find a similar double peak structure 
reveals that the internal velocity dispersions have a strong 
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Figure 5. Le/i: Evolution of the mean velocity dispersions for the two progenitors within the listed radii scaled by the values associated 
with the first ranked progenitor (<ri) at z = 0.2. The interpretation of the lines and line-styles is the same as in Fig. [4] Right: Evolution 
of the mean velocity dispersions for the two progenitors within the listed radii based only on the particles which were associated with 
the related progenitor at z = 1.2. The dispersions are scaled by the mean values associated with first ranked progenitor at z = 0.2. 
The interpretation of the lines and line-styles is the same as in Fig. [J] 



impact on the total kinetic energy. The structural similarity 
of the initial peaks in the left and the right panel of Fig. 
indicates that the increase in velocity dispersion is mainly 
caused by an internal acceleration due to the deepening po- 
tential well during peri-centre passages. This statement is 
supported particularly by the evolution of the velocity dis- 
persions within the 25/i _1 kpc spheres. As shown in Fig. [3] 
the distance of the cores is ~ 30/i _1 kpc resulting in only a 
partial overlap of the spheres associated with the two pro- 
genitors. Therefore, the enhancement of the total velocity 
dispersion is based mainly on the original particles of the 
associated progenitor and should resemble the behaviour in 
the right hand panel, as it does. However, a slight differ- 
ence between the evolution of the total velocity dispersions 
and those based only on original progenitor particles can be 
seen for almost all initial peaks. The peak associated with 
the second ranked progenitor (thin lines in the right panel) 
is slightly higher than the peak for the particles descending 
from the first ranked progenitor. For the central spheres this 
trend is more pronounced during the subsequent peri-centre 
passages, resulting in substantially larger velocity disper- 
sions for the particles descending from second ranked pro- 



genitor after total coalescence (thin lines on the right panel 
of Fig. [5). 

Combining the results derived for the evolution of the 
densities with the findings of this section results in the fol- 
lowing picture. After completion of the merging process the 
densities based on the particles descending from the sec- 
ond ranked progenitor are lowered relative to the densities 
related to the first ranked progenitor. The velocity disper- 
sions behave in the opposite way, they are higher for the 
particles descending from the second ranked progenitor. The 
progenitors have nearly equal masses but different concen- 
trations, which motivates the conclusion that an equal mass 
merger transfers a larger fraction of the released gravita- 
tional energy to the system of lower concentration and/or 
the energy is distributed amongst fewer particles in the low 
concentrated sy stem. Th i s beh aviour is probably related to 
the findings by iDehnenl (120051 ) who showed based on the 
phase space analysis of merging systems that the steepest 
cusps survive. However, it is not clear whether there exists a 
corr elation between the c entral slope and concentration (see 
e.g. iDiemand et~aT]|2004 ). 
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Figure 6. Left: Evolution of the mean velocity anisotropy parameter j3 for the two progenitors within the listed radii. The velocities of 
all the particles within the corresponding radii are taken into account. Heavy lines display the evolution of the first ranked progenitor and 
thin lines represent the second ranked progenitor. Right: Evolution of the mean velocity anisotropy parameter /3 for the two progenitors 
within the listed radii. Only the velocities of the particles associated with the related progenitor are considered. Heavy lines display the 
evolution of the first ranked progenitor and thin lines represent the second ranked progenitor. 



3.3 Velocity anisotropy parameter 

If during the relaxation process the particles of the lower 
concentrated system are pushed further away from the 
central regions this may have an impact on the veloc- 
ity anisotropy of the particles descending from this second 
ranked progenitor. The left panel of Fig. [6] displays the evolu- 
tion of the velocity anisotropy parameter (3 = 1 — 0.5 a t /oy , 
where ay indicates the tangential velocity dispersion and oy 
is the radial velocity dispersion, for all particles within the 
corresponding spheres. At z = 0.7, the time of the first peri- 
centre passage of the progenitor cores, the velocities within 
500 and 1000 h~ kpc appear to be momentary isotropic 
which is caused by the large tangential contribution of the 
passing cores and does not reflect the overall behaviour of 
the velocity fields. For smaller inclusion radii the anisotropy 
drops during and a short time after the second peri-centre 
passages. After z = 0.4 a stable configuration is achieved 
with a slightly radial anisotropy and a tendency to more 
isotopic velocities towards the centre. 

The velocity anisotropy based only on the particles de- 
scending from one of the progenitors at z = 1.2 (right panel 



of Fig. [5} support the above discussion of the velocity dis- 
persions. After the final coalescence of the progenitor cores 
the particles related to the second ranked progenitor show 
higher radial anisotropies than particles associated with the 
first ranked progenitor. The analysis of the densities and 
the velocity dispersions have already revealed that the par- 
ticles descended from the less concentrated progenitor are 
more strongly dispersed. To make that happen these parti- 
cles have to be provided with some additional radial velocity 
component which is supplied by the gravitational energy re- 
leased by of the merging systems. The result of this process 
is an relative increase of the radial anisotropy causing higher 
values of the /3 parameter. 

3.4 Entropy 

We define the entro py of the DM similar to the entropy of 
an id e al gas (see e.g. | White fc Naravan|[T987l ; iNavarro et al.l 
1 19951 ; lEke et al.l 119981 ). In analogy to the Sackur-Tetrode 
equation we write for the specific entropy of the dark mat- 
ter Sdm = ln(a^^/ pom) , where o"dm is the DM velocity 
dispersion and an additive constant is neglected. Follow- 
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Figure 7. Left. Evolution of the mean entropies {K = er 2 /p 2 / 3 ) for the two progenitors within the listed radii scaled by the values 
associated with the first ranked progenitor (K\) at z = 0.2. The interpretation of the lines and line-styles is the same as in Fig. [4] Right. 
Evolution of the mean entropy (K = <r 2 /p 2 / 3 ) for the two progenitors within the listed radii based only on the particles which were 
associated with the related progenitor at z = 1.2. The mean entropies arc scaled by the values associated with the first ranked progenitor 
(Ki) at z = 0.2. The interpretation of the lines and line-styles is the same as in Fig. [J] 



ing the convention used in the studies of galaxy clusters 
we d efine Kum = ""dmPdm' 3 ( see a l so iFaltenbacher et all 
120061 ). The DM specific entropy is clos ely related to the 
phase space density Qdm = Pdmo"^ (|Tavlor fc Navarrol 
l200ll ; lAscasibar et alj|2004h . 

Fig. [7] displays the mean entropies within the listed 
radii along the tracks of the two most massive progenitors. 
During the merging event the entropy within all radii rises 
temporarily. But it is difficult to estimate the change from 
initial to final entropy since the initial values are ill deter- 
mined, e.g. the central entropy of the first ranked progenitor 
(thick lines) shows a maximum value at 2 ~ 1 and drops by 
~ 20 per cent before dramatic growth takes place during the 
first peri-centre passage. In Fig. [S] we show evolution within 
the density-velocity dispersion plane between z = 0.85 and 
z = 0.2. The open circles indicate the starting point of the 
tracks and the filled circles give the location at z = 0.2. At- 
tached to both of these points is a dotted line indicating the 
loci of constant entropy or equivalently phase-space density. 
The distance between a pair of dotted lines gives the ratios 
between the final and initial entropies Kf/Ki, which range 
from 0.95 to 1.15 for the lOO/T^kpc and the 1000/i _1 kpc, 



respectively. This changes are rather small, in particular if 
one considers the occasionally sever deviations from isen- 
tropic behaviour during the active merging phase. Within 
the 500/i _1 kpc sphere initial and final entropies are equal 
(heavy long dashed line in Fig. [7}. This radius may be of 
particular interest, since it is the largest radius included by 
the virial radius for all times considered here. 

For the 100 and 250 fe~ 1 kpc spheres the evolution of the 
entropies show a clear double-peaked structure at the first 
peri-centre passage. This feature is most prominent for the 
100 /i _1 kpc sphere. The analysis of the the densities (Fig. UJ 
revealed that within this sphere the density during the first 
peri-centre passage increases strongest compared to the be- 
haviour of other spheres. This particular steep density en- 
hancement can counterbalance the increase of the entropy 
caused by the peaking velocity dispersion. 

We have shown that particles descended from the sec- 
ond ranked progenitor build a subset of particles within the 
final cluster with lower central densities and higher velocity 
dispersions compared to particles descended from the first 
ranked, more concentrated progenitor. The combination of 
these two processes causes a substantial rise in the central 
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Figure 8. Evolution of the first ranked progenitor for the listed inclusion radii in the density-velocity dispersion plane. The slanted 
lines show loci of constant entropy K = a 2 /p 2 ^ 3 of equivalently phase-space density p/cr 3 . The tracks begin at a redshift of z = 0.85 
(approximately the beginning of the merging event) marked by open circles and end at a redshift of z = 0.2 after the progenitor cores 
have merged indicated by filled circles. Units are arbitrary and relative values of the entropy increase toward the upper left. Compare 
with Fig. 8 of Navarro et al. (1995). 



entropies for the particles descended from the second ranked 
progenitor as displayed in Fig. [7] This is a further indication 
that the central parts of the less concentrated progenitor get 
heavily disrupted whereas the more concentrated progenitor 
is less affected by the merging event. 



3.5 Central Potential 

We compute the contribution of the various spheres to the 
central potential by direct summation. 



(<r) = -G J2 



JV(<r) 



r, 



(1) 



Here, r cen is the centre of the progenitor, m p and r% are the 
mass and the position of the ith particle within the inclusion 
radius r and G is the gravitational constant. The evolution 
for the central potential of the most massive progenitor is 
displayed in Fig. [5] In order to see how the evolution of the 
potential is affected by occasional very central passages of 
single particles we have computed the potential of the par- 
ticle which closest to the centre for the given moment. We 
find that the contribution of the most central particle to the 
potential is on average —0.003 x 10 6 ( kms -1 ) 2 with maximal 
amplitudes of ~ — 0.01 x 10 6 ( kms -1 ) 2 . Therefore, occasional 
central passages of single particles do not seriously affect the 
overall evolution of the gravitational potential. The oscilla- 
tory behaviour seen in the evolution of the density and the 
velocity dispersion is also apparent in the central poten- 
tial. Since the potential based on the 1000 /i _1 kpc results in 
the most pronounced variations we compare in Fig. [TU] the 
central potentials along the tracks of the two progenitors 
based on the 1000/i _1 kpc sphere. Obviously the central po- 
tential of the second ranked progenitor undergoes stronger 
variations during the oscillatory relaxation phase. The six 
peri-centre passages based on the core distances (Fig. [3} are 
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Figure 9. Impact of the merger on the central gravitational po- 
tential computed within the listed radii along the track of the 
most massive progenitor. 
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Figure 10. Evolution of the central gravitational potential com- 
puted within the 1 /i _1 Mpc centred along the tracks of the two 
most massive progenitors. The solid line displays the evolution 
of the first ranked, more concentrated progenitor and the dotted 
line shows the behaviour of the second ranked, less concentrated 
progenitor. 

clearly visible in the evolution of the central potential. The 
time interval of the later oscillations is ~ 10 8 yr. 



4 SUMMARY AND DISCUSSION 

We have analysed the merging of a cluster-sized dark mat- 
ter haloe focusing on the evolution of the density and the 
velocity dispersion within fixed proper radii. This approach 
is complementary to the more common analysis based on 
the scaling relations for dark matter haloes. Whereas the 
latter emphasises the self similar nature of gravitationally 
collapsed systems, the former provides insight into the pro- 
cesses at fixed scales and may help to understand the cause 
of the self similarity. Our main results are: (1) About 3 Gyrs 
elapse from the first encounter of the progenitor outskirts 
to the final coalescence of their central cores. Within the 
~ 2 Gyrs between the first peri-centre passage and the final 
core merging six peri-centre passages take place. During that 
process the periods decrease by a factor of ten from ~ 10 9 yrs 
to ~ 10 8 yrs. (2) The central densities of the merging pro- 
genitors peak at each peri-centre passage and are most pro- 
nounced at the first one. The density enhancement is caused 
by the overlapping volumes of the subhaloes but also by the 
contraction due to the deepening of the potential well dur- 
ing a peri-centre passage. (3) After the merger is completed 
(z < 0.4) the densities and the velocity dispersions within 
the fixed proper radii keep roughly constant. The increase 
for larger radii at z < 0.1 is associated with an ongoing 
intermediate merger as di scussed in § O In agreement with 
iRomano-Diaz et al.l (|2006h we find that the properties of the 
dark matter haloe change rapidly during merging events but 
remain constant during phases of passive evolution. (4) The 
ranking of the progenitor concentrations determines the fi- 
nal distribution of the particles. Less concentrated haloes get 
more dispersed and show higher radially anisotropic velocity 
dispersions. 

Finally, we briefly discuss the impacts of the oscilla- 
tory relaxation process of the collisionless matter on the 
baryonic component. It is of interest to consider the pos- 
sibility that enhanced AGN activity in the cluster-centred 



galaxies may be stimulated by the large oscillations of the 
haloe potentials during a major cluster merger. We note that 
the period of relaxation oscillations for our merging cluster, 
~ 10 s yrs, is similar t o the duty cycle observed in AGNs 
jVoit fc Donahue! 120051 ). although only a small fraction of 
AGNs reside in massive clusters. Nevertheless, to explore 
this idea, we performed a simple gasdynamical calculation 
to determine how much additional radiative cooling occurs 
when the merging potentials are deepened, compressing the 
hot cluster gas. For this purpose we set up a traditional 
quasi-steady state cluster cooling flow in a potential of mass 
3 x 10 14 with typical concentration c = 6.74. The gas tem- 
perature at large radii, 3.5 keV, w as chosen according t o the 
observed cluster M — T relation (|Arnaud et all l2005h and 
the initial (hydrostatic) gas density was adjusted until the 
baryon mass matched the cosmic value at the virial radius. 
The relaxed steady state co oling flow that results c losely re- 
sembles the observations of IVikhlinin et all (|2006l ) and has 
a central cooling rate of M = 43 Mq yr -1 (assuming abun- 
dance 0.4 solar). The potential was then set into sinusoidal 
oscillations with period 0.6 Gyr and an amplitude similar to 
that of Figure 10 by sinusoidally varying the concentration 
between c and 0.43c. During these oscillations, the cooling 
rate M(t) near the centre of the flow varied with an am- 
plitude of only 7 per cent, much less than the unperturbed 
steady state value which is at least ~ 10 times larg er than al- 
lowed by X-ray observations (|Peterson et aljfeOQl ). We con- 
clude that the additional radiative cooling stimulated by the 
oscillating merger potentials is negligible. 

Perhaps a more fruitful approach along these lines is 
the possibility that galaxies (and large black holes) orbit- 
ing about the cluster core can be processed by dynamical 
friction into the central galaxy/black hole during the relax- 
ation process. The acceleration due to dynamical friction 
(du/dt)df oc mp/u 2 depends on the mass m and the veloc- 
ity u of the orbiting object and the density of the ambient 
material. By setting the velocity equal to the mean disper- 
sion u = er, we see that the time scale for dynamical fric- 
tion tdf = u/(du/dt)df oc cr 3 /p = A' 2//3 . During apo-centre 
passages when the dark halo entropy passes a minimum, ob- 
jects orbiting near the cluster-centred galaxy may efficiently 
reduce their angular momentum, increasing the chance of a 
capture by the central black hole during the subsequent con- 
traction phase. It is conceivable that such events, stimulated 
by the global merger, could result in enhanced AGN activity 
with a period imprinted by the dark matter oscillations. 
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